(N 
> 
O 

oo 
(N 

(N 
O 



>< 



ON DEGENERATE PLANAR HOPE 



o 
o 

(N 

Q : BIFURCATIONS 



^ : M. R. Ricard 

Q 



Faculty of Mathematics and Computer Science, 
University of Havana, 

CP. 10400, Cuba 
rricard@matcom.uh.cu 

Abstract 



H ' Our concern is the study of degenerate Hopf bifurcation of smooth 

planar dynamical systems near isolated singular points. To do so, 
we propose to split up the definition of degeneracy into two types. 
Degeneracy of first kind shall means that no limit cycle surrounding 
the steady state can emerge after or before the critical point, with 
the possible emergence of limit cycles surrounding the point at infin- 



ity. Degeneracy of second kind shall means that either several limit 
cycles or semistable cycles as a limiting case, emerge surrounding the 
steady state super or sub critically. In degenerate bifurcation of second 
kind we also show that the radius of the emerging cycle tends to zero 
with an "anomalous" order as the bifurcation parameter tends to the 
critical value. Finally, we give a sufficient condition for degenerate 
bifurcations of second kind up to 6-jet-equivalence, and show some 
"typical" forms for degenerate bifurcations. 

AMS Subject Classification: 34C23, 37G15, 34C29 
Keywords: Hopf bifurcation; limit cycles; averaging method 

1 Introduction 

The goal of this paper is the study of degenerate Hopf bifurcation (HB) near 
the critical value of a one-parameter family of smooth planar vector fields in 
a neighborhood of an isolated singular point. The HB is usually associated 
to the emergence of hmit cycles after or before the critical value of the bifur- 
cation parameter. An isolated periodic orbit of a vector field on the plane, is 
a limit cycle. The HB phenomenon is studied recurrently in literature, and 
its main details and implications for dynamical systems of higher dimensions 
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via the Center Manifold Theorem are gathered in the monograph [T] . Several 
efforts to compute generalized HB in finite dimensional systems can be found 
in literature (see [2] and references therein). Applications of the HB in finite 
dimensional dynamical systems often arise in the study of population-based 
models, for instance in predator-prey systems featured by the Allee effect [3]; 
in chemical systems, as in the Schnakenberg's [1]; in the study of coupled 
systems near a supercritical HB [5J; in the theory of electronic circuits [6]; 
in mathematical economics [7j; in the modeling of mechanical systems [8]; 
even in the study of travelling waves phenomena [H] and, generally speak- 
ing, almost whenever the considered model shows non-linear oscillations. 
Remarkable applications of HB in fluid dynamics can be found in [1] and 
references therein. In addition, the topic is extensively treated in the theory 
of bifurcations with higher codimension, as in the Turing-Hopf instabilities 
|10j . Nevertheless, in order to focalize the scenario, we suggest in this paper 
to split the classification of degeneracy into two different kinds. Degeneracy 
of first kind shall means that no limit cycle surrounding the steady state 
emerges neither after nor before the critical value of the bifurcation param- 
eter, with the possible emergence of limit cycles surrounding the point at 
infinity [11]. The second kind of degeneracy shall means that multiple limit 
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cycles or semistable cycles as a limiting case, emerge either at super- or at 
subcritical bifurcation surrounding the steady state. We focus our attention 
in such a system that, in a neighborhood of the origin, can be written in the 
form 

X = F{X,a) (1) 



where 



R 



I . \ 



k+l=l 



a 



is the Taylor expansion of F in a neighborhood of the origin, the coefficients 
cr^ = cr^ (a) are real, and X = {x,y)'^ . Here a is a real parameter, and F 
is at least a C'^-function of a. Any smooth dynamical system in a vicinity of 
an isolated singular point Pa on the plane, can be rewritten in the form Eq{T] 
by a parameter-dependent shift of coordinates which relocates the singular 
point Pa of the system into the origin of coordinates. Eq{T] is called the 
system in variations corresponding to a given former system near the singular 
point Pa- Any reference to the steady state Pa of the former system will be 
associated to the punctual orbit at the origin of the Eq{TJ The vector field 
with polynomial components of degree R represented in the main part in 
Eqj2l is called the R-jet of EqO] around the origin, which is usually denoted 



+ (2) 



Jr (F) (0). Two smooth vector fields defined in a neighborhood of the origin 
are called R-jet equivalent if their i?-jet coincide. 

Dynamical systems featured by an HB often arise in mathematical models 
for biological or chemical reaction systems [1]. In general, these models lead 
to polynomial vector fields having fifth degree at most. Moreover, some 
complex interactions between the different components in the system may 
lead to a representation in which a rational fraction with linear or quadratic 
denominator appear in the former system. For instance, this is the case in the 
Michaelis-Menten kinetics or in the prey-predator Rolling systems in presence 
of AUee effect [3]. In such situations, after rescaling the temporal variable, 
the system can be transformed into an orbitally equivalent [12] polynomial 
system of degree at most six. 

We recall that only few combinations of coefficients corresponding to odd 
terms in the reaction system can contribute to the emergence of limit cycles 
at bifurcation. We have called these the Hopf coefficients. These Hopf co- 
efficients can be easily calculated and will play a similar role than such the 
Lyapunov coefficients do [T3]. We suggest a new classification of degenerate 
bifurcations in terms of the Hopf coefficients, through the "discriminant" in- 
troduced in [To]. In this direction, we introduce the concepts of degenerate 
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HB of first and second kind. 

In accordance with the proposed classification we subdivide the HB Theo- 
rem, in order to emphasize differences in the resulting behavior at bifurcation. 
We prove that, at degenerate bifurcation of first kind no limit cycle surround- 
ing the steady state emerges, being allowed the emergence of limit cycles 
surrounding the point at infinity. At supercritical (subcritical) degenerate 
HB of second kind multiple or semistable limit cycles emerge surrounding 
the steady state. In correspondence to the procedure in [10], we shall see 
that, at a non-degenerate HB, a single limit cycle emerges. Each kind of 
degenerate HB is complemented with appropriate examples. 

The period of the emerging limit cycles in one-parameter bifurcations at- 
tract the attention of researchers (see [TT]) showing that, the main term in 
the asymptotic of the period of the emerging periodic solution, characterizes 
the bifurcation. Here we focus our attention in a procedure that gives us 
simultaneously an asymptotic estimate of the radius of the cycle and the 
frequency of the corresponding periodic orbit. As we shall see, the radius of 
the limit cycles and the period of the periodic solutions emerging at degen- 
erate HB of second kind, tends to zero with an "anomalous" order as the 
bifurcation parameter tends to the critical value. 
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The plan of the paper is as follows. In Section [2] we summarize previous 
results in the treatment of an HB taken basically from |10j. In Section [3l 
the notion of discriminant is quoted, which is used in Definition [3] to classify 
HB taking into account its asymptotic behavior as the trace of the Jacobian 
tends to zero. The splitted version of the HB Theorem in [H [2] and [3] in this 
Section describe the emergence of limit cycles in non-degenerate, degenerate 
of first and of second kind respectively. Finally, in Section H] are gathered 
"typical" forms corresponding to each type of degenerate bifurcation. We 
include a sufficient condition for degenerate HB of second kind of Eq{T] valid 
up to 6-jet-equivalence . 

2 Preliminaries 

Let 5a and be the determinant and the trace 



5, 



a 



det (J,) 



(3) 



a 



trace {J a) 



(4) 
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of the Jacobian matrix 




J, 



a 



(5) 




at the origin O of the function F in Eq{T] The subindex in Eqsj3] and H] 
indicates a functional dependence respect to a, which varies in an open small 
neighborhood U of the point a* at which the trace vanish to change its sign. 
We shall assume that, the function a is an homeomorphism between U 
and a neighborhood of = 0, so the transversality condition {a^) 7^ 
is not required. Consequently, the parameter Tq, \Ta\ ^ 1, can be considered 
as the intrinsic bifurcation parameter (see [TU]) and, = is the critical 
value. The HB appears provided the inequality 



holds for any value Ta in some neighborhood of Tq = 0. Then, for every value 
a G t/, or equivalently, for every value Ta G V , we have 



rl - 45„ < 



(6) 



5a >0. 



(7) 
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Further, if < (respect, > 0) the origin is a stable (respect, unstable) 
focus. The bifurcation is subcritical if there is a limit cycle emerging for 
negative values of close enough to zero. The bifurcation is supercritical 
if there is a limit cycle emerging for positive small values of r^. From EqE] 
follows (Tq^-ct^q < and we do not loose generality assuming (Jqi < 0. Besides, 
there is no added restriction if we consider the Jacobian matrix in the system 
Eq{T] to have the simplest form: 



Ja 2 



1 



(8) 



where 

A, = +v/4(5a - r2 > 

in a neighborhood of Ta = 0. So, we shall assume in the following that 
the Jacobian in EqlS] has the form EqlHl If necessary, the system Eq|T] can 
be rewritten through a linear transformation of variables in order that the 
Jacobian has the required form Eq|8], whenever the condition Eqj6] holds. 
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2.1 Averaging Hopf periodic solutions 

Let us quote in this Subsection some results about the procedure in the study 
of HB proposed in [TU]. Following that paper, we rewrite the system EqlH 
in the form 

X = J„X + ^ (X) (9) 

where X = {x (t) ,y (t))"^, and the vector function \E' is given by the Taylor 
expansion of the difference 

^ (X) = F (X, a) - JaX 

so, (X) contains all nonlinearities. In the function ^ (X) we include the 
Taylor terms until a required precision, say up to i?-jet-equivalence, together 
with the corresponding remainder. First we assume that (X) is analytical 
and a bit later, in Remark [1], we turn back into the smooth case. In [TO] , 
the authors proposed an algorithm allowing the reduction of an analytical 
reaction system showing an HB into a second order differential equation 
representing a weakly nonlinear oscillator in normal form. This transform 
of variables is analytical and nonlinear in general, but it was proved there 
that it is enough to consider the linear part of this transform to obtain the 
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equation of the oscillator preserving the required accuracy. The main idea of 
this procedure is that the transform of variables can be taken "close" to the 
appropriate linear transform in a neighborhood of the origin. So we quote the 
procedure in [TU], and consider an invertible analytical transform of variables 
between neighborhoods of the origin 

Y = n{x) = TX + g (x) (10) 

where T = {'Jij) is a non-singular matrix, and Q (X) be analytical without 
linear terms. By the Inverse Function Theorem, the existence of the in- 
verse is guaranteed because F is non-singular and Ti has smooth continuous 
derivatives. The inverse to Eq{TO]has the form 

X = (Y) = T-^Y + IC{Y) . (11) 

Definition 1 We say that the diffeomorphisms Ti and F have a contact at 
the origin of order 5 G N, 5 > 1, 

n{x)-vx = o{\\xf^^'^ (12) 
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as \\X\\ 0. 

We would like to substitute the function H in Eq{TO] by an equivalent 
simpler one, say F, in the procedure of rewriting the system Eq{T] in new 
easy-handled variables to be used in Proposition [H In this endeavor we get 
more precision as the order of contact S between 7i and F be greater. Further, 
from EqIT2] follows that exists whenever F is invertible, and it can be also 
concluded that H^^ and F^^have a contact of order S. Following the ideas in 
[To] it can be proved that, if F satisfies certain "concordance" condition and 
7i has a contact with F at the origin of order R, then Y = Ti. {X) represents 
an analytical transform of coordinates such that every solution (x (t) ,y (t)) 
to the analytical system Eq{T] is transformed into the form 



being z (t) unknown. Hence, the integration of the system Eqj9] can be re- 
duced to the integration of a second order differential equation in the variable 
z. Consequently, rather than the exact expressions of the functions Q {X) 
and }C{Y), it is sufficient to take the linear transforms F and F~^ instead 
of Ti. and Ti.^^ in the derivation of the equation in z, considering that this 




(13) 
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substitution preserves the required accuracy. More precisely, for a given R 
in Eql2] we may take a transform of coordinates Ti. which has a contact of 
order R with F at the origin. In the following statement we quote a result in 
Pm Prop.l] in which are gathered all the above ideas: 

Proposition 1 Let us assume that Eql^i holds. Then, there exists an in- 
vertible analytical transform of variables EqlT^ in the system Eql^ such that 
Eq[T^ holds. The matrix T is any non-trivial linear combination of the pair 



/ \ 

1 

.1 ^1 






1 








^01 



(14) 



The function z in EqlWi satisfies the following second order equation: 



Z — TaZ + 5a Z = G {^Z, Z 



(15) 



where the right hand side in Eq\T^ does not involve linear terms in z, z. More 
precisely, 

G {z, z^ = n2 {r ( j,/c {Y) + ^ {n-^Y)) (le) 
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being II2 the standard projector over the second component. 
So, the function in Eq{T6] can be expanded asymptotically by 

G'(^,^) =n2{r(^(r-ir))}+o(||Ff+^) . (i7) 

Let us now return to EqJTSl We shall look for an oscillation with positive 
and small, but finite, amplitude e. The small parameter e is connected with 
the small bifurcation parameter and will be defined a bit later. Taking in 
EqJTSl the change of variables 

z (t) = {t) (18) 

we get the equation of a weakly nonlinear oscillator in normal form: 

- + 5a = e G (^q, . (19) 

Then, to each periodic solution to Eql9]will correspond a non-trivial periodic 
solution to EqUS In [TD], was considered the Krylov-Bogoliubov averaging 
method [H] to derive an asymptotic expansion to the solution of EqJTOl To 
do so, let us consider the new variables r = r (t) and 9 = 9 {t) defined as 
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follows 

— r cos {t + 9) 
q — —rsm{t + 9) 

then, the corresponding averaged equations are 

2tt 



1 f 

r = — — / sin0 {— Ta r sin 4> + sG {r cos 4>,—r sin 4>;e)} d4> 

2'K J 







27r 



9 — —- — / COS0 {— Ta r sin (p + eG (r cos 4>, —r sincp] e)} 
ZTrr J 



thus, 



where 



r ^^{ra-p{r;e)} 



9 ^q{r;e) 



2-K 

p(r;s) = — / sin 4> G (r cos 4>, —r sin 4>; s) d(f) 
7rr J 



27r 

qir]e) — —- — / cos ^ G (r cos 0, — r sin 0; s) d(f) 
ZTrr J 
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Let us quote now some important properties about the functions p and q 
above. 

Proposition 2 Functions p (r; e) and q (r; e) have at least order O (e^) or, 

equivalently, p{r;e) /r^ and q{r;e) /r^ have a finite limit as r ^ 0. More- 
over, the Taylor expansions ofp (r; e) and q (r; e) must not contain odd powers 
ofr. 

From Proposition [2] it can be concluded that the development of p has 
the form 

p (r; e) = ps e^r"^ + p^ e^r^ H (28) 

in which ps = Ps (Ta), as a consequence of the dependence that the right 
hand side of Eq{T]has on the bifurcation parameter Ta- In Subsections 14. 
14.21 and 14.31 can be found some examples showing this type of dependence. 

We recall now the fact that, in accordance with EqsI^H and [2HI the M 
coefficients {p2j+i)jLi wiU determine the bifurcation up to 2 (M + l)-jet- 
equivalence. So, it can be expected that we would manage the coefficients 
P2j+i by taking M appropriate independent relations involving the parame- 
ters in the system EqlH This scenario is called a codimension-M bifurcation. 
In accordance with our purpose in this paper, let us now introduce a defini- 
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tion that will be essential for the classification of HB. The reason to take it, 
will arise in the next Section. 

Remark 1 We also recall that we are considering smooth vector fields, while 
in the construction of the analytical transform of coordinates in EqlWi we 
use analytical properties. But notice that the transform EqlT^ as well as the 
appropriate R-jet approximation, which is analytical, lead to functions EqlT7\ 
in which the principal part remains unaltered after their substitutions into 
the smooth system EqUl 

Definition 2 Let p2s+i be a coefficient in the formal development Eq^El 
which is derived from the formal oo-jet of F. It shall be called negligible if 
satisfies 

\p2s+l\<Ks\Ta\ (29) 

for certain constant Kg > as Ta —>■ 0. The function p{r\e) in Eq\2^ is said 
to be negligible if for all s & N the coefficient p2s+i is negligible. 

Naturally, P2s+i = and p = are respectively included in the above 
definition. For instance, p = if in the formal joo {F) (0) the non-vanishing 
terms have even degree. Moreover, we can get P2s+i = in spite of the 
existence of non-zero coefficients with degree 2s + 1 in the formal Taylor 
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development of the right hand in Eq|T] As we shall see in the next Section, 
negligible terms have no influence in the generation of limit cycles. Besides, 
with this notion we are able to give a more detailed version of the Proposition 
3 in PO] as follows. 

Proposition 3 // the function p (r; e) is non-negligible, there must exist a 
positive integer N and a positive real value = Tq (tq) such that p (r, e) has 
the non-trivial Taylor expansion: 



where a; = +1 or —1. In addition, the behavior of the factor as Ta ^ 
obeys the following alternative: either 



'^2N+2 ^2N+2 



(30) 



lim r. 



-2N 




= r 



2N 



> 



(31) 



or, for a given 7, < 7 < 1 




(32) 
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As in [12], the symbol O5 in EqI32] means a sharp estimate, that is: 
= O {\Ta\^) and r^^^ ^ o{\Ta\"') as ^ 0. We recall that the bifur- 
cation is supercritical (respect, sub critical) if a; = +1 (respect, a; = — 1). In 
the supercritical case, the root tq appears for Ta > (respect, < 0), so 
the limit in EqlST] or the order relation in Eq. [32] should be considered as 
Ta — > 0+ (respect, — > 0—). Consequently, 

Proposition 4 Let us assume that p (r; e) be non-negligible and also, that 
ro in Eql3^ has the property in Eq[31\ Hence, there is a positive root p to the 
equation 

p (r; e)-Ta = (33) 

either for positive or for negative values Ta close enough to zero. Further- 
more, up to the leading term, the root to Eq\3^ has the form 

P=\^) (^* + 0(|r.|)) +0(6^) . (34) 

Proof. The property EqJ3T]is equivalent to Tq"^^ = r~^^ + (|Ta|) and, also 
to ro = r, + 0(|r„|). ■ 

Due to Proposition H] the function p in Eql26] has been called a discrim- 
inant for the HB in [TU]. Let us assume the existence of a (finite) positive 
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root EqlM] such that EqlHI] holds. Then, the small parameter e introduced 
in Eq{T8] can be taken as 

= \Ta\ . (35) 

From Eql35] and the relation \Ta\ = O ^|Ta|^^^j as Ta ^ for > 1, it 
follows that EqlSH can now be written as 

P = n+ 0(|rj^/^) . (36) 

If EqJST] does not hold we may assume Eql32| so in a similar way as we do 
in Proposition H] to get EqJ53l we now arrive to 

e'"" = |r.r^ . (37) 



Moreover, if for instance 



'^O = '^L kaT + O (kaD aS 



for certain positive number tl, then EqJHUcan be rewritten as 



I |(l-7)/2Af I^L , ^ / 2\ 

P=\ra\ ^Oie) 
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and follows 

P = r,+ O (|rj(^--)/^) . (38) 

The hypothesis of the HB Theorem, as appear in [TOl Theorem 1], con- 
tains an implicit reference to the non-degenerate case, more precisely, to the 
condition EqjSTl As we shall see in Section [3l EqsjSTl [32] and [29] put in 
evidence the reason for differentiation in HB we are suggesting here. In ac- 
cordance with the intention and terminology in this paper, we shall give in 
Section [3] a split version of the HB theorem. 



3 Theorems for degenerate Hopf bifurcation 



Basically, the standard classification of HB is conformed by two main classes: 
while the system at degenerate HB shows a center at the critical value a^., 
the system at a non-degenerate HB shows a weak focus at the critical value, 
leading to the emergence (super- or subcritical) of a limit cycle [3]. When 
the bifurcation occurs in a one-parameter family of vector fields whose first 
non-zero derivatives at the origin have order 2N + 1, > 1, it is called 
[TT] a generalization of the Andronov-Hopf 's. Other higher codimension HB, 
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for instance the Bautin bifurcation, are often called generalized [12]. In 
our formulation, we shall include such higher codimension HB, hence we 
shall assume we have a family of systems parametrized by the bifurcation 
parameter Tq. Consequently, the coefficients in Eq|2H] also depend on the 
bifurcation parameter, ps = Ps (Ta), so we can a priori classify the HB taking 
into account how these dependences are. 

Definition 3 We shall say that the Hopf bifurcation is degenerate of first 
kind if p is negligible (See Definition\E )■ Let N be given as in Eq\3^ The 
bifurcation shall be called degenerate of second kind, if there exists a number 
7, < 7 < 1, such that EqlSEholds. The HB shall be called non-degenerate, 
provided Eq[31[ 

For instance, if p = in Eq l25l the bifurcation shall be degenerate of first 
kind. Further, the existence of at least one non-negligible P2s+i derived from 
the formal cxD-jet of F, no matter how large the number s is, implies that the 
HB will not be degenerate of first kind. We also remark that, at degenerate 
bifurcation, we are implicitly considering that the system moves close to a 
higher codimension point as the bifurcation parameter varies, because the 
non-zero coefficients p2s+ivanish at = 0. 
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Remark 2 There is a reason to take < 7 < 1 m the definition of the 
degeneracy of second kind above. As F in EqUl is a continuous function of 
Ta in a vicinity of Ta = 0, we may assume •y > as a consequence of Eq[3E. 
Furthermore, from Eq^^ and the "smallness" of the parameter e in Eq\T^ as 
Ta — ^ follows that 1 — 7 > 0. 

It is easy to see, from Definition [3l that degenerate (of any kind) HB 
implicitly implies that the system shows a center at the critical value of the 
bifurcation parameter. This is because the main terms of p in Eqj2l] tends 
to zero as — > 0, so the origin is a center at the critical value. We split up 
the HB theorem taking into account each type of bifurcation. We recall that 
the bifurcation is supercritical if the cycle emerges provided < <C 1, and 
subcritical if it emerges for < — ^ 1. 

Theorem 1 (non-degenerate Hopf bifurcation) Let us assume that Eq\^ 
holds and that Eq\3^ has a root Eq\3^ with the property Eql3J\ for positive 
(respectively, negative) but sufficiently close to zero values of the bifurcation 
parameter Ta- Then, a single limit cycle to the system in Eql^ emerges. 
Furthermore, the limit cycle is orbitally asymptotically stable (respect., un- 
stable) if and only if the bifurcation is supercritical (respect., subcritical). 
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The radius of the emerging cycle is r = Os [\Ta\ j, while the frequency is 

ZU = 1 + O ^|Ta|^^^j as Ta — 0. 

Proof. The proof follows from Eql2S] and Eql^H Let the main term in 
Eqj30] have the property Eqj3Tl Then, a single root to Eqj33] tends to zero 
as Tq — > 0. This fact means that a single limit cycle emerges at bifurcation. 
It is not excluded the existence of further different roots to Eql33], but if any 
other appears, it determines a limit cycle that does not vanish in a vicinity of 
Ta = 0, SO the cycle "persists" along the bifurcation. As follows from EqsJMj 
|55] and [56l we get the order of the radius of the limit cycle. The order of the 
frequency is determined in Subsection 13. 3[ ■ 

Consequently, a necessary but not sufficient condition for the emergence 
of a limit cycle at bifurcation, is the existence of nonzero odd order terms in 
the expansion of the components in the right-hand side of EqUJ For instance, 
the following result states that the condition is not sufficient. 

Theorem 2 (degenerate HB of first kind) Let us assume that EqlB\holds 
and assume p negligible (Definition\^. Then, neither limit cycle surrounding 
the steady state emerges after nor before the critical value. 
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Proof. From Eq|2n] and Eq|2S] follows that none of the roots to Eq|33] tends 
to zero as 0. More precisely, if any root exists, it tends to infinity. 

Taking into account only the leading terms we may write 



P2k+i = ak rl+''^ 



where /x^ > 0, are constants which can be zero. Let corresponds to the 
first non-zero a^, and take (if the bifurcation is supercritical), so 

EqM is 



CXD OO 



Mfe+{n-iV)/iV^27V 



n=N n=N 

hence, if there is a root to Ta — p (r; £:) = then it should tend to infinity as 



Remark 3 We remark that, at a degenerate HB of first kind, further limit 
cycles may persist in a neighborhood of the critical value of the bifurcation 
parameter as can be seen in Subsection \4-3\ Furthermore, if the order of a 
negligible term p2s+i is greater than one, then limit cycles surrounding the 
point at infinity may emerge in the so-called HB at the infinity (see [Tl^ 
and references therein). In Subsection \4-3\ we shall show an example of this 
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situation. For instance, this is the case if 



P2S+1 




) 



being Ta 0, (3 > 0, A ^ 0. 

Now, with the same procedure in the proof of Theorem [H we have, 

Theorem 3 (degenerate HB of second kind) Let us assume that Eqldi 
holds and that EqlS^ has a root EqlMi with the property Eq\3^ for positive 
(respectively, negative) but sufficiently close to zero values of the parameter 
Ta- Then, it can be assured the emergence of at least one limit cycle to the 
system in Eq\^ the radius of which has order 



< 7 < 1 corresponds to the one in EqUE. 

Proof. The proof for the non-degenerate bifurcation can be repeated, but 
taking the small parameter from Eqj371 Notice that in this case more than 
a single positive root to EqI2S] tending to zero as — > 0, may appear. ■ 




(39) 




The number 7, 
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At a degenerate HB of second kind, different behaviors can be observed: 
a single or several limit cycles may emerge, including semi-stable cycles as 
the limiting case at which different limit cycles collapse. In Section H] some 
examples with this kind of degeneracy are considered to show that, multiple 
limit cycles or semistable limit cycles might appear either sub- or supercrit- 
ically. Bearing in mind the construction of the examples in Section HJ looks 
easy to find sufficient conditions for the emergence of multiple limit cycles, 
as it is done in Subsection 14.51 

We can show sufficient conditions for the stability of the emerging limit 
cycles surrounding the steady state. Such conditions are based on the be- 
havior of the discriminant p near the root to which the cycle corresponds. 

Proposition 5 Consider a root p to Eq [^ so it corresponds to a limit cycle 
L. This cycle is asymptotically stable or unstable if the number dp/dr (p) is 
negative or positive respectively. 

Proof. In the calculation of dp/dr (p) we can assume p = ro as in Eq|3Sl 
or p = r^, in Eqj38l in accordance with the type of bifurcation. Due to 
the continuity argument near = in Eqs. or [3H1 and the fact that 
dp/dr (ro) (or dp/dr {ri)) does not vanish, the assertion follows. ■ 
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3.1 The Hopf coefficients 



Let us take M > 2 in Eq{T]and the map Eq{TO]in which Ti has a contact with 
r of order M at the origin, we get 



l<k+l<M 



a 



\Y' 



where T = {'Jij) and T ^ = (jJ^ij). For instance, if we take T = Fi, the 
quantities 

/^ii = 1 ; ^12 = ; f^2i = ^aK^ ; ^^22 = -2A^^ (40) 

are the components of F^^. So, EqJT6] can be written 



G{z,z 



^ Rki {^^^llZ + /iia^;) (/^2i^ + ^^22z) (41) 

2<A;+KA/ 



+o((||r| 



A/+1 



where 



= 721 C^L + 722 



(42) 



Hence, Eqj4T] can be rewritten in terms of powers of the F-components 
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as 



g(z,z)= J2 H,,z'' (z) +0((\\Y\f^')) . (43) 



2<k+l<M 

Only the non-zero Hmn in EqH3] corresponding to such pairs (m, n) for which 
Kmn ^ 0, where 

/•27r 

Kmn = / COS™ sin"+^ # (44) 
Jo 

can contribute to the appearance of a nonzero term in Eql30] so, to the 
appearance of a hmit cycle solution. 
Let us introduce the following 

Definition 4 We shall call the Hopf coefficient of degree {2N + 1) to the 

coefficient p2N+i in the expansion Eq WBf which is an algebraic combination 
of coefficients Hmn in Eq\4^ provided m + n = 2N + 1. 



From EqMl EqgS] and EqI2S] follows directly that 



Ps = (3iJo3 + (45) 



is the Hopf coefficient of third degree, and 



p, = -l{5Ho5 + H23 + H^i) (46) 
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is the Hopf coefficient of degree five. Other Hopf coefficients can be derived 
by forward calculations. We recall that are required only the coefficients of 
r and the coefficients of the main part of the Taylor expansion, to calculate 
the Hopf coefficients. Becomes easy to check from EqSl] (see Appendix E]) 
that the coefficients in EqJl3] leading to Hopf coefficients in Eq{T7]up to 6-jet 
equivalence are 

(47) 

-f^0,5 -f^2,3 -f^4,l ■ 

In accordance with EqH2| if the Jacobian matrix have the form Eq|8] and 
being F = Fi, we have the numbers 

Rmn = \ {ra t^^,„ " K (^In^n) (48) 

We recall that the matrix F, as a linear combination of the matrixes in EqJT^ 
depends on the elements of the matrix which is supposed to have the form 
in EqJHl Then, taking F = Fi and F^^ = (/ij^) , we get the following equalities 

i703 = -(8i?03)Aa' , (49) 
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H21 = -2 (i?21 A' + 2i?12 TaAa + 3i?03 T^) A^^ , (50) 
def 

Ho, = - (32 i?05) (A.-') , (51) 

H23 = -8 (i?23 A^ + 4i?i4 TaA, + 10i?05 T^A^) A'^ , (52) 
def 

H^, = -2 {R^, At + 2Rs2 rAi + 3i?23 r^A^ (53) 

def 

+ARu rlA^ + 5i?o5 r^) A;^ , 

from EqH3l 

It is concluded in this Subsection that, up to 6-jet-equivalence, we can 
easily check the type of degeneracy the degenerate HB shows. Particularly, 
EqsJlS] and |l6] give the full information for polynomial vector fields of degree 
at most six. In Section H] we show different examples of degenerate HB within 
the class of polynomial vector fields of degree at most six. 

3.2 Asymptotic expansions of the limit cycles 

Going back to the substitutions given in EqUH] and EqUUl we could derive 
uniform asymptotic expansion to the solution of Eq{T] So, if the HB is non- 
degenerate, it is possible to develop the periodic solution (t) = {x (t) ,y{t)) 
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to Eq{T] generating the limit cycle, by 



X{t) = Ui{t) i\Ta\)^^ + O {\Ta 



(54) 



y{t)=viit) (|r„|)2^ +Oi\Ta 



(55) 



and 



with frequency 



Ul{t) 

it) 



r* r 



-1 



COS [w t) 
^ — sin {tu t) ^ 



w = l + q(p,\Ta 



(56) 



(57) 



and period T = From Eq|25] the angular speed of the oscillation is 
obtained. Note that, up to the leading terms, the expansions of the cycle 
solution in EqsJS^ and [55] are uniform, as the O ^ | | ^ -terms are bounded 
functions. The expansions in EqsHH and [55l can be taken also for the limit 
cycle that emerges at degenerate HB of second kind considered in Theorem 



[31 but substituting the exponent ^ by and r^, by the given in Eqj38l 
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3.3 On the period of the hmit cycles 

Here we give a brief comment about the period of the hmit cycles, considering 
the interest in the topic [TT]. The period of the emerging hmit cycles in 
non-degenerate or degenerate of second kind HB can be determined from 
the formula for the frequency, given in EqlSTl To do so, it is necessary to 
consider EqsJ2^ [271 and Up to 6-jet-equivalence, the expansion results 

q = —e^p' (^[3/730 + H,^] + \e^p' [5/^50 + ^32 + ^u]) + O [e^) . 

The coefficients H^i in the above formula, can be calculated in a similar 
way as the others in EqsHU to [53l If the HB is non-degenerate we have 
the relation Eql35l while we have to consider Eq|37] if the bifurcation is 
degenerate of second kind. Further, the factor p take the values Eql3S] or 
Eqj38] in accordance with the type of bifurcation. The procedure yields to 
the following estimates of the frequency as — > 0: 

9 = C'(|rJ^) (58) 
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for the emerging cycle in the non-degenerate case, and 

g = C'(|r,|^) (59) 

if the bifurcation is degenerate of second kind. We remark the fact that, in 
the last two formulas, the order it is not necessarily sharp. In the examples 
in the next Section we also have done a reference about the period. 

4 Typical forms in degenerate HB 

The main concern in this Section are degenerate HB of first and second kinds. 
In Subsections l4.1l to l4.4l we shall study "typical" forms by using polar coordi- 
nates in order to give a familiar, more geometrical, description of the nature 
of the bifurcation. Notice that the examples in the referenced Subsections are 
not normal forms, because in these cases no genericity conditions (see [12] ) 
are involved. It is not difficult to show that, using the averaging method 
we can get the same conclusions. In Subsection 14.51 we state, as an example 
of what can be expected in presence of a degenerate HB of second kind, a 
sufficient condition using the Hopf coefficients up to 6-jet-equivalence. 
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4.1 Multiple cycles in supercritical degenerate Hopf 

bifurcation of second kind 

In this paragraph we shall give an example of a system which shows a de- 
generate HB of second kind in accordance with Definition [S] 



X 



-1 



v 



1 



J 



( \ 



( \ 



( \ 



2\2 



X 



(60) 



Neither limit cycle surrounds the focus while the trace r = 2a^ is negative. 
At the critical value a = 0, the origin is clearly a center to the system Eqj60l 
Further, when the trace becomes positive the steady state becomes unstable 
and two small limit cycles emerge with radius 



respectively. Both cycles emerge due to the bifurcation and have radius 
Os (t^^^) as r ^ 0. The first cycle is stable, while the second is unstable. 
To check the above assertions we only need to rewrite the system EqjGni in 
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polar coordinates. We get the system: 



I 



r = ar 



From this example it can be concluded that, at degenerate bifurcation of 
second kind, several limit cycles may emerge. Of course, with the same idea, 
it is possible to build polynomial dynamical systems with higher degree 2A^+1 
showing the emergence of different limit cycles at a super- or subcritical 
degenerate HB of second kind. The frequency of the emerging cycles is 



as iV = 1 and 7 = 2/3 in EqEH 

4.2 Semistable cycles in supercritical degenerate Hopf 

bifurcation of second kind 

The following system shows a supercritical degenerate bifurcation of second 
kind: 
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/ . \ / , \ / \ 

-1 



X 



\ 



1 



J 



X 



-2c? [x^ + y^) 



( \ 



( \ 



2\2 



+a [x^ + y^) 



X 



(61) 



The origin is an stable focus without any surrounding hmit cycle while the 
trace r = 2a^ is negative. At the critical value a = 0, the origin is clearly 
a center to the system EqEH The system shows a degenerate bifurcation 
of second kind in accordance with Definition [31 Further, when the trace 
becomes positive the steady state becomes unstable and, a single small limit 



cycle emerges with radius 



r = a 



1/2 



This limit cycle is the cu-limit set of any orbit inside the circle with the 
exception of the origin, but it is the a-limit set of any orbit outside the 
circle. The corresponding system in polar coordinates is: 

r = ar (r^ — a)^ 



9 = 1 



37 



From this example it can be concluded that, at a degenerate bifurcation, 
semistable limit cycles may emerge. The radius of the limit cycle is r = 
Os {t^^^) as r ^ 0. The frequency of the emerging cycle is 

as = 1 and 7 = 2/3 in Eqj59l 



4.3 Degenerate Hopf bifurcation of first kind without 



emergence of limit cycle 



1. Let us first consider a system showing a degenerate HB of first kind: 



X 



a —1 
1 a 



( \ 



X 



( \ 



X 



(62) 



The origin is an stable focus without limit cycle surrounding while 
a < 0. At the critical value a = the origin is a center, and the steady 
state turns unstable if a > 0. It can be noted the existence of a limit 
cycle with radius r = 1, but this cycle is not a consequence of the 
bifurcation because it persist along the bifurcation. The polar system 
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is: 




showing that the cycle r = 1 is a stable limit cycle for a > 0, changing 
its stability in dependence of the sign of a. This circle is still an orbit 
of Eqilfor a = 0. 

2. We recall in the fact that, the persisting limit cycle may be semistable, 



as in the system 



X 



a —1 
1 a 



-2a {x^ + y^) 



X 



( \ 



X 



\' J 
(63) 



which also shows a degenerate bifurcation of first kind. The corre- 
sponding polar system is: 



r = ar (1 — r 



9= 1 



SO, the limit cycle r = 1 is semistable, its interior stability changes with 
the sign of a, and it is still an orbit of EqESlfor a = 0. 



39 



4.4 Degenerate HB of first kind showing limit cycles 
at infinity 

Let us consider now a system showing a degenerate HB of first kind, but 
leading in this case to the so called HB at infinity [TT]: 



/ . \ / 

X 



a —1 
1 a 



\ 



( \ 



X 



( \ 



X 



(64) 



P > 0. Here we take = sign (a) ■ \af . The portrait becomes clear taking 
polar coordinates: 

r = ar (l — r^) 

e = i 

We conclude that, as a consequence of the degenerate bifurcation of the 
system Eqj6l] no limit cycle emerges surrounding the origin, but an unstable 
limit cycle surrounding the point at infinity, whose radius is r = a~^/^, 
emerges supercritically. More precisely, the cycle emerges for small positive 



values of the trace t = 2a. 
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4.5 A sufficient condition for degenerate Hopf bifurca- 
tion of second kind 

We shall consider here a degenerate HB for the system Eq{T] up to 6-jet- 
equivalence. The following statement, which is inspired in the examples 
in Subsections 14.11 and \A.2\ gives a sufficient condition for the existence of 
degenerate bifurcations of second kind. The idea in this assertion is very 
simple. 

Proposition 6 Let Eql^i holds for the system EqUl at the origin, and let 
Jq (F) (0) be such that the coefficients ps in Eq^^ and in Eq^^ satisfy the 
conditions: 

P3 = r:Qs + o {t2) ; p5 = -rl^-' Q, + o (r^^"^) 

for some 'j (1/2 < 'j < 1) as Ta —>■ 0, where Q^, are both positive numbers 
and 

A = Ql~4Q5>0 . 

Then, two different limit cycles if A > 0, or one semistable limit cycle if 
A = 0, emerge at the supercritical degenerate bifurcation. The radius of the 
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cycles are Os [r^a ^^^"^j as Ta — * 0+. 

Proof. Taking 1/2 < 7 < 1 in Eq|32] follows EqJSTl hence e = ri^^^^^^ Up 
to the O (tq) leading term in Eql3Sl we get 

Q5 - Qs + 1 = 

meaning that the algebraic equation above has either two different positive 
roots {k = 1, 2), or a single positive root with multiplicity two, depending 
on A. Then, the proof follows from Eqj33]and EqlMl ■ 

5 Conclusions 

To study degenerate Hopf bifurcation in smooth dynamical systems near an 
isolated singular point on the plane we first focus in a classification of such 
bifurcations via a discriminant function. If the bifurcation is non-degenerate, 
then a single limit cycle emerges. The radius of this cycle can not be a pri- 
ori supposed to have order O {ja^'^^ as — > 0, but in general have order 
O ^r^^^ j , for some integer A^. In this scenario, the emerging limit cycle may 
coexist with other cycles which persist along the bifurcation. This classifica- 
tion gives some light to the question of whether degenerate bifurcations lead 
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to the emergence of limit cycles. We show that, no limit cycle surrounding 
the steady state emerges neither at super- nor at subcritical HB when it is 
degenerate of first kind. In this scenario, limit cycles surrounding the point at 
infinity may also emerge. For a degenerate HB of second kind, we found the 
anomalous asymptotic order of the radius and period of the emerging limit 
cycles, and further, we give a sufficient condition to the appearance either 
of a couple of limit cycles or one semistable cycle of the system Eq{Tl up to 
6-jet-equivalence. Finally, we propose some examples, which are called here 
"typical" forms, showing different behaviors that can occur at supercritical 
degenerate HB, being either of first or of the second kind. 
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6 Appendix 

In this Appendix we include the values of the integrals in EqjHl which are 
used in the implementation of the Krylov-Bogoliubov averaging method. Fur- 
ther, we include the coefficients Hmn leading to Hopf coefficients. Notice that, 
for m + n <Q the non-zero numbers are: 



For a given matrix F, which is a non-trivial linear combination of the pair 



K21 = in ; 






in Eq{Tll we denote F ^ 



{fj'ij)- Then, the Hmn leading to Hopf coefficients 



up to 6-jet-equivalence are 



21 



3-^30/^11/^12 + 2-R2l/^ll/^12/*21 + -^21/^11/^22 



+2-Rl2/^ll/^2l/*22 + -^12/^12/^21 + 3-Ro3/^2l/*22 ' 



Ho3 — -^12/^12/^22 + -^03/^22 + -^30/^12 + -^21/^12/^22 
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Hii — 2i?32/ii]^/i2l/^22 + 5-R50Atii/^l2 + -^41/^11/^22 + 4-R4l/ill/^12/^21 
+2i?23/^ll/^12/*21 + 3-R23/ill/^2l/^22 + 3-^32/^11/^12/^21 
+-^14/^12/^21 + 4-Rl4/^ll/^2l/^22 + 5-Ro5/^2l/^22 ) 

2 3 3 2 2 

-f^23 = 10i?50/iii/ii2 + 4-R4i/iii/ii2/i21 + 6-R4l/^ll/^12/^22 

2 3 2 2 2 

+6-R32/iii/ii2/i2l/^22 + -^32/^12/^21 + 3-R32/ill/^l2/^22 

2 2 3 2 2 

+6i?23/^ll/^12/^2l/i22 + -^23/^11/^22 + 3-R23/^12/i2l/^22 

2 2 3 2 3 

+6-Rl4/^12A''2l/^22 + 4-Rl4/ill/i2l/^22 + 10-^05/^21 /^22 ) 

-f^05 = -^32/^12/^22 + -^23/^12/^22 + -^14/^12/^22 
+i?05/^22 + R50IA2 + -^41/^12/^22 ) 

where i?^; are given in EqH2l 
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